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ABSTRACT 

Trypanosoma brucei causes human African 
trypanosomiasis and regularly switches its major 
surface antigen variant surface glycoprotein (VSG) 
to evade mammalian host immune responses at 
the bloodstream form (BF) stage. Monoallelic 
expression of BF Expression Site (BES)-linked 
VSGs and silencing of metacyclic VSGs (mVSGs) in 
BF cells are essential for antigenic variation, 
whereas silencing of both BES-linked and ml/SGs 
in the procyclic form (PF) cells is important for cell 
survival in the midgut of its insect vector. We have 
previously shown that silencing BES-linked VSGs in 
BF cells depends on T6RAP1. We now show that 
Tib RAP 1 silences both BES-linked and ml/SGs at 
both BF and PF stages. The strength of T6RAP1- 
mediated BES-linked VSG silencing is stronger in 
the PF cells than that in BF cells. In addition, 
Formaldehyde-Assisted Isolation of Regulatory 
Elements analysis and MNase digestion dem- 
onstrated that depletion of TbRAPI in PF cells led 
to a chromatin structure change, which is signifi- 
cantly stronger at the subtelomeric VSG loci than 
at chromosome internal loci. On the contrary, no 
significant chromatin structure changes were 
detected on depletion of TbRAPI in BF cells. Our 
observations indicate that TbRAPI helps to deter- 
mine the chromatin structure at the insect stage, 
which likely contributes to its strong silencing 
effect on VSGs. 

INTRODUCTION 

Human African trypanosomiasis is caused by infection of 
Trypanosoma brucei and is inevitably fatal without treat- 
ment. Inside the human host, bloodstream form (BF) 



T. brucei cells stay in extracellular spaces and express 
variant surface glycoprotein (VSG) as its major surface 
antigen that is exposed to the host immune system (1,2). 
To evade the host's immune responses, T. brucei cells 
undergo antigenic variation and regularly switch their 
VSG coat (3), which is essential for a persistent infection. 
Although there are >1000 VSG genes and pseudogenes in 
the T. brucei genome (4), in BF cells, VSGs are expressed 
exclusively from BF VSG Expression Sites (BESs), which 
are polycistronically transcribed by RNA Pol I (5). The 
VSG gene is the last gene in the BES and is located 
adjacent to the telomere while the BES promoter is 
usually 40-60 kb upstream (6). The T. brucei 427 strain 
used in this study has ~20 nearly identical BESs (7), 14 
of which carry distinctive VSG genes (8). However, at any 
moment, only one BES promoter is fully active, resulting 
in a single type of VSG being expressed. VSG monoallelic 
expression ensures the effectiveness of antigenic variation 
and is essential for T. brucei virulence. 

Trypanosoma brucei is transmitted through its insect 
vector, tsetse (Glossina spp.). After ingestion by tsetse 
into its midgut, T. brucei differentiates into the procyclic 
form (PF) and expresses procyclins as surface glycopro- 
teins, of which the C-terminus is resistant to protease 
cleavage (9). VSG is susceptible to protease degradation, 
and all VSG genes are silent in PF cells. Therefore, VSG 
silencing at the insect stage is critical for T. brucei to 
survive in its environment, but the underlying mechanism 
is poorly understood. Once T. brucei cells migrate into the 
salivary glands of tsetse, they differentiate into the 
metacyclic form, shed procyclins, express metacyclic 
VSGs (mVSGs) on the cell surface and acquire infectivity 
again. Each cell expresses one type of mVSG, but the 
population expresses many different mVSGs (10). This 
heterogeneity will presumably facilitate the population 
to establish an infection in mammalian hosts (11). 
Similar to BF VSGs, mVSGs are also transcribed from 
subtelomeres by RNA Pol I, but transcription of mVSGs 
is monocistronic (12,13). mVSG expression is silenced 
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within a few days after the initial infection, allowing sub- 
sequent antigenic variation to take place, but the mechan- 
isms of mVSG silencing have not been extensively studied. 

Accumulating data have revealed that multiple mechan- 
isms are involved in the regulation of BES-linked VSG 
monoallelic expression at the BF stage, including limited 
accessibility to RNA Pol I (14), restricted transcription 
elongation (15,16), chromatin remodeling (17-22) and 
telomeric silencing (23). Telomeres, located at the ends 
of linear chromosomes, are essential for genome stability. 
In several organisms including yeast, human, mouse, 
T. brucei and Plasmodium falciparum, telomeres form a 
heterochromatic structure that suppresses the transcrip- 
tion of sub telomeric genes (24-29). Trypanosoma brucei 
telomere DNA consists of thousands of duplex TTAGG 
G repeats (30). We have previously identified 77?TRF as 
the duplex telomere DNA-binding factor (31) and 
77? RAP 1 as a 77?TRF-interacting partner (23). Depletion 
of 72? RAP 1 led to derepression of all subtelomeric silent 
BES-linked VSGs and simultaneous expression of multiple 
VSG proteins on the cell surface at the BF stage (23). 
However, whether 72? RAP 1 has a similar VSG silencing 
function in PF cells and the mechanism of 72? RAP 1- 
mediated silencing are not clear. 

In this study, we found that depletion of 72? RAP 1 by 
RNA interference (RNAi) led to derepression of mVSGs 
at both the BF and PF stages. Silencing of BES-linked 
VSGs in PF cells depends on 72? RAP 1 and is stronger 
than that in BF cells. Furthermore, removal of 72? RAP 1 
led to a loosened chromatin structure in PF cells, particu- 
larly at the subtelomeric VSG loci, but not significantly in 
BF cells, indicating that 72? RAP 1 -mediated silencing 
involves modulation of chromatin structure in PF cells. 

MATERIALS AND METHODS 

Reverse transcription and quantitative reverse 
transcriptase-PCR 

Total RNA was extracted using RNAstat 60 (Tel-Test 
Inc), treated with DNase I (Qiagen) and purified using 
RNeasy columns (Qiagen). Reverse transcription, 
quantitative reverse transcriptase-PCR (qRT-PCR) and 
quantification were performed as in (23). 

Formaldehyde-assisted isolation of regulatory elements 

Formaldehyde-assisted isolation of regulatory elements 
(FAIRE) analyses were performed according to (32). 
Cells were fixed with formaldehyde before the chromatin 
was sonicated in a Bioruptor (Diagenode). Free DNA was 
extracted by phenol chloroform. Input samples were 
reverse cross-linked before DNA was extracted. The 
amount of DNA was estimated by quantitative PCR 
using iTaq SYBR Green Supermix with ROX in an 
Opticon II (Bio-Rad), and the amount of FAIRE-extracted 
DNA was normalized with that of the input DNA. 

Micrococcal nuclease digestion 

In all, 5 x 10 7 cells per sample were permeabilized with 
digitonin and treated with 1 unit of Micrococcal 



nuclease (MNase) (Worthington Biochemicals) for 1, 2.5 
or 5 min. DNA was isolated from treated cells, separated 
by agarose gel electrophoresis, blotted onto a nylon mem- 
brane and hybridized with specific probes. 

Statistical analysis 

P-values of unpaired Mests were calculated using 
GraphPad Prism. 

RESULTS 

7#RAP1 is essential for normal cell growth in PF cells 

To examine the functions of 77? RAP 1 in PF cells, we 
introduced an inducible 77? RAP 1 RNAi construct in 
29-13 cells that expresses the Tet repressor and the T7 
polymerase (33). Several independent clones were 
obtained that exhibited different growth defects on induc- 
tion of 77?RAP1 RNAi (Supplementary Figure S1A-C). 
To reduce the phenotypic variations often associated with 
clonal cell lines, we obtained a pool of 77? RAP 1 RNAi 
cells (referred to as 'PRi-pool') for subsequent detailed 
analyses, and critical phenotypes were confirmed in the 
PRi-C2 clone (Supplementary Figures SIB, D and F and 
S8B). 

On induction of 77? RAP 1 RNAi, we observed a 
decreased 77? RAP 1 protein level and a severe growth 
defect by day 2 in PRi-pool cells (Figure 1A and B, 
Supplementary Figure S1E), indicating that 77? RAP 1 is 
essential for cell viability. Similar 77? RAP 1 knockdown 
and growth arrest were seen in PRi-C2 cells (Sup- 
plementary Figure SIB and F). The control 29-13 cells 
carrying an empty RNAi vector did not show any 
growth defects or decrease in the 77? RAP 1 protein level 
when treated with the same amount of doxycycline 
(Figure IB, Supplementary Figure S1E). Flow cytometry 
analysis was performed to examine cell cycle profile after 
depletion of 77? RAP 1. We observed a mild decrease in the 
Gl population and a mild increase in the G2/M and the 
sub Gl populations in PRi-pool cells (Figure ID, 
Supplementary Figure S2A), which is similar to what we 
observed in 77? RAP 1 -depleted BF cells in our current and 
earlier studies (Supplementary Figure S2B and C) (23). 

BES-linked VSG silencing depends on TbRAPl at the 
insect stage 

We compared the steady-state mRNA levels of several 
BES-linked VSGs before and after depletion of 77?RAP1 
by qRT-PCR. Thirty-five percent of 77?RAP1 was still 
present at day 1 after induction of RNAi, but <10% of 
77? RAP 1 was left by day 2 (Figure 1A). We therefore 
examined VSG derepression at day 2 and day 3. All 
tested VSG mRNA levels increased several 10- to several 
100-fold after depletion of 77?RAP1 (Figure 2 A, see 
Supplementary Table SI for VSG nomenclature). 
Western analysis also detected a prominent amount of 
VSG2 in PRi-pool cells (Figure 1A, Supplementary 
Figure S1E). Therefore, 77? RAP 1 is essential for silencing 
BES-linked VSGs in PF cells. Similar to that in BF cells, 
expression of RNA Pol II-transcribed 77?PGI (a glycolytic 
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Figure 1. Depletion of TbRAPl led to a cell growth defect and derepression of BES-linked VSG2 in PF cells. (A) Western analysis of various protein 
levels in PRi-pool (left) and T^TRF RNAi (right) cells. Total cell lysates were prepared at days 0, 1, 2 and 3 after adding doxycycline. (B) Growth 
curves for PRi-pool and control cells with the empty vector under induced (+ Dox) or un-induced (—Dox) conditions were plotted as Population 
Doublings versus days after induction. Average Population Doubling values were calculated from six independent experiments. Error bars represent 
standard deviation. (C) Gating schematic of un-induced PRi-pool cells in the flow cytometry analysis. Sub-Gl cells contain DNA contents less than 
2C and often experience DNA degradations. PI, propidium iodide. (D) Quantification of populations of PRi-pool cells at different cell cycle stages 
before (day 0) and after (days 1, 2 and 3) depletion of TbRAPl. Unpaired t- tests were done to compare day 1, 2 and 3 values with the day 0 value. 
Significant differences are marked with asterisks. *, 0.01 <.P<0.05; **, 0.001 <P<0.01. 



protein), 7&RPS15 (a ribosomal protein) or RNA Pol 
I-transcribed rRNA was not affected, indicating that 
72>RAPl's effect on VSGs is specific. In addition to 
BES-linked and mVSGs, VSG genes located on minichro- 
mosomes (which mainly consist of internal 1 77 bp repeats 
and terminal telomere repeats) are also at subtelomeres 
but lack upstream promoters. Silencing of a 
minichromosomal VSG, VSG671, was not affected by de- 
pletion of TbRAPl — we did not detect its mRNA in 
northern blotting (data not shown), although the 
VSG671 gene was detectable in Southern analysis 
(Supplementary Figure S3) — indicating that VSG dere- 
pression requires a functional upstream promoter. 
Similar VSG derepression phenotype was also observed 
in PRi-C2 cells (Supplementary Figure SID). On the 
contrary, no VSG derepression was detected in control 
cells by qRT-PCR (Supplementary Figure S4A) or 
western blotting (Supplementary Figure S1E). In 
addition, induction of TbTKF RNAi led to a cell growth 
arrest within 2 days in PF cells (31), but not VSG2 expres- 
sion (Figure 1 A), indicating that VSG derepression is not a 
consequence of cell growth arrest. Although 7&RAP1 
interacts with TbTKF (23), it has not been proved that 
TbTRF is strictly required for localizing 7&RAP1 at telo- 
meres. Therefore, depletion of TbTRF does not necessar- 
ily have the same VSG derepression phenotype. 

J^RAPl is required for silencing of m VSGs in BF and 
PF T. bvucei cells 

mVSGs are normally silent in BF and PF cells, but the 
silencing mechanisms are poorly understood. In our 



laboratory strain 427, five mVSG genes have been 
identified: mVSGs 397, 531, 639, 653 and 1954 (34). 
qRT-PCR analysis showed that mVSG mRNA levels 
increased several to several 10-fold on depletion of 
72? RAP 1 (Figure 2B), whereas no significant changes 
were detected in control cells (Supplementary Figure 
S4B). We performed similar analyses in BF cells using 
two cell lines carrying pZJM(3-based (35) 72? RAP 1 
RNAi vectors. One includes the 72? RAP 1 BRCT 
fragment (Ri-BRCT) and the other the TbRAPl C- 
Terminus (Ri-CT). Both cell lines exhibited growth 
arrest and derepression of BES-linked VSGs on induction 
of 77? RAP 1 RNAi (Supplementary Figure S5), indicating 
that these cells behaved similarly to the Ri-2 and Ri-9 
76RAP1 RNAi cells studied previously (23). mVSGs 
were derepressed several to over a 1000-fold on 72? RAP 1 
depletion in Ri-BRCT and Ri-CT cells, but not in control 
cells (Figure 2C; Supplementary Figure S4C and D). 
Therefore, Tb RAP 1 is required to silence mVSGs in 
both BF and PF cells. 

Comparing the J^RAPl-mediated VSG silencing at 
different life cycle stages 

We found that 76RAP1 silences BES-linked and 
metacyclic VSGs, but the derepression of BES-linked 
VSGs appeared to be much stronger in PF cells than in 
BF cells. In our previous study, on Tb RAP 1 depletion, 
VSG 3, 6 and 11 were derepressed up to 8-, 23- and 41- 
fold, respectively, in BF cells (23), whereas the same genes 
were derepressed up to 270-, 220- and 260-fold, respect- 
ively, in PF cells (Figure 2A). 
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Figure 2. Depletion of 77? RAP 1 resulted in derepression of BES-linked VSGs in PF cells (A) and derepression of mVSGs in PF (B) and BF (C) cells. 
Steady-state mRNA levels for BES-linked VSGs 2, 3, 6, 11, mVSGs 397, 531, 639, 653 and control genes at days 0, 2 and 3 (PF) or at 0, 24 and 36 h 
(BF) after adding doxycycline in T&RAP1 RNAi cells were estimated using qRT-PCR and normalized against that before induction. Day 0 value is 
set to 1 but not shown. The fold change in mRNA levels for three independent induction experiments was shown. 
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Figure 3. Tb RAP 1 -mediated BES-linked VSG silencing is stronger in 
PF cells than in BF cells. The relative mRNA levels for the same VSG 
genes were estimated by qRT-PCR when the VSG is in BF or PF 
TbRAPl RNAi cells before (silent) and after (derepressed) 7#RAP1 
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Ri-2 cells (expresses VSG2) (23), and VSG2 derepression was analyzed 
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Single Marker (SM) expresses VSG2 (33), pVS3-2/ODl-l expresses 
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To better quantify VSG derepression, we used 
qRT-PCR to carefully compare individual BES-linked 
VSG mRNA levels at the silent state in wild-type (WT) 
PF and BF cells, at the derepressed state in 7&RAP1- 
depleted PF and BF cells, and at the fully active state in 
BF cells. The steady state VSG mRNA levels under 
various conditions were normalized against that in PF 
WT cells, which was arbitrarily set to 1. Three BES- 
linked VSGs — VSGs 2, 3 and 9 — were analyzed owing to 
availability of all necessary cell lines in our laboratory 
(Supplementary Table S2). We found that first, in WT 
cells, VSGs 2, 3 and 9 mRNA levels were 28-, 16- and 
9-fold higher in BF cells than in PF cells, respectively 
(Figure 3, PF WT and BF WT). This is likely owing to 
the fact that the VSG mRNAs are less stable in PF cells 
than in BF cells (36). Second, 7&RAP1 was essential for 
VSG silencing at both stages. Depletion of 7&RAP1 
elevated the VSG2, 3 and 9 mRNA levels on average 
333-, 527- and 217-fold, respectively, in PF cells and on 
average 22-, 133- and 30-fold, respectively, in BF cells 
(Figure 3, PF / 76 RAP 1 RNAi -/+ Dox, BF/ 
7Z>RAPlRNAi -/+ Dox). Third, the derepression of 
BES-linked VSGs resulting from 76 RAP 1 depletion was 
significantly stronger in PF cells (217-527-fold) than that 
in BF cells (22-133-fold). The P- values were all <0.002 for 



VSG2, 3 and 9, in unpaired Mests when comparing 
changes in BF and PF cells. Therefore, the 76RAP1- 
mediated BES VSG silencing was stronger in PF cells 
than in BF cells. We noticed that in un-induced 
PRi-pool cells, the VSG mRNA level is -3-fold lower 
than that in PF WT cells (Figure 3), possibly reflecting 
small variations between different cell populations. 

In yeast, longer telomeres confer stronger telomeric 
silencing (37). However, telomeres in our PF and BF 
76 RAP 1 RNAi cells are of similar length (on average 
~15kb), suggesting that telomere length is not the 
reason for different Tb RAP 1 -mediated silencing in BF 
and PF cells. To rule out the possibility that 72? RAP 1 
might not be equally depleted in BF and PF cells, which 
could contribute to the differential FSG-derepression at 
the two stages, we examined 72? RAP 1 protein levels in 
all 76RAP1 RNAi cells. In Ri-2 and Ri-9 cells, 76RAP1 
was depleted as efficiently as described previously (23). At 
the time points when various 76 RAP 1 RNAi cells were 
examined, comparable residual amounts of 76 RAP 1 were 
detected (on average, 4% in Ri-2, 14%— 1 1 % in Ri-9 and 
16% in PRi-pool; Supplementary Figure S6A-C) (23). 
Therefore, it is unlikely that different levels of VSG dere- 
pression in BF and PF cells were due to different levels of 
76RAP1 depletion. Another possibility is that 76RAP1 
was expressed at different levels in BF and PF 
cells. However, by northern and western analyses, we 
found that 76 RAP 1 mRNA and protein levels were ap- 
proximately the same at the two life cycle stages 
(Supplementary Figure S6D and E). 

Although 76 RAP 1 is also required for mVSG silencing, 
a higher degree of variation in mVSG derepression was 
observed in independent inductions of Tb RAP 1 RNAi 
(Figure 2B and C), and no significant differences were 
observed between BF and PF cells (P > 0.05), except 
that derepression of mVSG397 in Ri-BRCT cells is 
slightly higher than that in PRi-pool cells (P = 0.03). 

Depletion of 7#RAP1 led to more loosely packed 
chromatin structures in PF cells 

To further investigate whether Tb RAP 1 -mediated VSG 
silencing involves a chromatin remodeling mechanism, 
we compared the subtelomeric chromatin structure 
before and after depletion of 76 RAP 1 using the FAIRE 
analysis (32), where DNA fragments free of bound 
proteins are enriched in the final phenol-chloroform- 
extracted fraction. We used FAIRE to extract 
subtelomeric BES-linked VSG DNAs, whose amount 
was estimated by quantitative PCR, and the fold change 
was calculated by dividing the post- Tb RAP 1 depletion 
amount by the pre-depletion amount. 

Depletion of 76 RAP 1 for 1.5 days did not cause a 
significant change in FAIRE values in PRi-pool cells 
(data not shown), probably owing to insufficient 
76RAP1 depletion (35% of 76RAP1 protein was left by 
day 1, Figure 1A). We therefore performed FAIRE at day 
2.5 in PRi-pool and found an 8-19-fold increase in BES- 
linked VSG DNA and a 4-7-fold increase in mVSG DNA 
enrichment, which were significantly higher than that in 
the control cells (Figure 4). To determine whether this 
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Figure 4. Depletion of 77? RAP 1 led to loosened chromatin structure, 
particularly at subtelomeric VSG loci, in PF cells. The amount of 
FAIRE-extracted DNA after depletion of TbRAPl for 2.5 days was 
quantified by PCR using primers specific to indicated gene loci and 
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was calculated from at least three independent experiments. Standard 
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Tb927.2.2440; Tbl510: Tb09.211.1510. *, 0.01 < P< 0.05; **, 
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effect is telomere-specific, we also examined chromatin 
structure at several random chromosome-internal gene 
loci, including 77>927.2.2430, 77>927.2.2440 and 
7*09.21 1.1510. Depletion of 77>RAP1 in PF cells led to 
a 3-4-fold increase in FAIRE-extracted DNAs at these 
loci, which was mildly but significantly higher than that 
in the control cells (Figure 4), suggesting that 77? RAP 1 
may have a more global effect on chromatin structure 
than we anticipated. However, the changes at 
subtelomeric VSG loci were significantly stronger than 
that at the control loci (P<0.05, Figure 4), indicating 
that the effect of 77? RAP 1 on chromatin structure is 
more prominent at subtelomeric loci. The BES promoter 
regions that are 40-60 kb upstream of the telomere were 
similarly affected by 77? RAP 1 depletion as the chromo- 
some internal regions (Figure 4). 

To rule out the possibility that the chromatin structure 
change resulted from elevated transcription rather than 
directly from loss of 77>RAP1, we examined the 
VSG671 locus, which is located at the subtelomere of a 
minichromosome (Supplementary Figure S3) and not 



derepressed on 77? RAP 1 depletion. Depletion of 
Tb RAP 1 led to an —5-fold increase in FAIRE extracted 
VSG671 DNA, which was significantly different from that 
in control cells (Figure 4), indicating that the VSG671 
chromatin structure was affected by 77?RAP1, and that 
this effect was not due to elevated transcription of 
VSG671. 

We subsequently performed the same FAIRE analyses 
in BF cells in multiple independent 77? RAP 1 RNAi cell 
lines (Ri-CT, Ri-BRCT, Ri-2 and Ri-9) at 24 and 36 h 
after 77? RAP 1 depletion. A mild increase in the amount 
of FAIRE-extracted BES-linked, metacyclic and 
minichromosome VSG DNA was visible in all 77? RAP 1 
RNAi and control cells at both time points (Sup- 
plementary Figure S7). A similar change was also 
observed for the chromosome internal genes and the 
BES promoter loci (1-3-fold) (Supplementary Figure 
S7). Unpaired Mest showed that all VSG loci behaved 
similarly to chromosome internal loci, and all 77? RAP 1 
RNAi cells behaved similarly to the control cells 
(P>0.05). Therefore, we did not detect significant 
changes in chromatin structure on depletion of 77? RAP 1 
at the infectious stage. 

To rule out the possibility that inefficient depletion of 
77? RAP 1 was the reason for the lack of significant chro- 
matin structure change in BF Tb RAP 1 RNAi cells, we 
examined 77? RAP 1 protein levels by western analysis. At 
24 h after 77?RAP1 RNAi induction, the 77?RAP1 protein 
level typically reduced to 6 or 26% of the WT level in Ri- 
CT and Ri-BRCT cells, respectively (Supplementary 
Figure S6C). By 36 h, Ri-2 and Ri-9 cells had -10% or 
less 77? RAP 1 left (Supplementary Figure S6A and B). 
Although in PRi-pool cells, Tb RAP 1 was typically 
depleted to —16% of the WT level after 2.5 days of induc- 
tion (Supplementary Figure S6C). Therefore, the different 
phenotypes we observed in BF and PF cells did not result 
from unequal levels of 77? RAP 1 depletion. 

To further validate our conclusion, we examined the 
accessibility of chromatin by MNase before and after de- 
pletion of 77? RAP 1. PRi-pool cells were induced for 2.5 
days and treated with a fixed amount of MNase for an 
increasing length of time. Subsequently, the genomic 
DNA was isolated and separated on agarose gel. 
Without MNase treatment, nearly all of the genomic 
DNA was in the limiting mobility size range (Figure 5, 
EtBr, dots). After 2.5-5 min of digestion with MNase, a 
portion of the chromatin was digested into mono-, di-, tri- 
nucleosomes and so forth, and the isolated DNA gave a 
laddering pattern with the smallest bands of —150, 350 
and 500 bp (Figure 5, EtBr, asterisks). After depletion of 
77?RAP1, the chromatin was more accessible to MNase 
digestion, and less material was left undigested while 
more DNA was in the mono-nucleosome fraction 
(Figure 5A, lanes 7 and 8; Figure 5B, lanes 9 and 10). 
The samples from 77? RAP 1 -depleted cells also appeared 
more smeared, whereas samples from un-induced cells 
gave sharper bands (Figure 5 A, compare lanes 4 and 8; 
Figure 5B, compare lanes 4 and 9), indicating that in 
77? RAP 1 -depleted cells, DNA was less protected by nu- 
cleosomes, resulting in more DNA degradation. Depletion 
of 77? RAP 1 altered the MNase digestion pattern mildly in 
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Figure 5. Depletion of TbRAFl led to chromatin structure changes in PF cells. Trypanosoma brucei cells were treated with MNase for an increasing 
length of time (min) before genomic DNA was isolated, separated on agarose gel, blotted onto nylon membranes and hybridized with specific probes. 
The PRi-pool cells were treated and hybridized with the VSG18 (A) or the TTAGGG repeat probes (B). The Ri-CT and Ri-BRCT cells were treated 
and hybridized with the VSG2 (C) or the TTAGGG repeat probes (D). Black dots and the arrow head represent limiting mobility in ethidium 
bromide (EtBr)-stained gels and hybridized blots, respectively. Asterisks and triangles represent mono-, di-, tri- and tetra- nucleosomes on EtBr- 
stained gels and hybridized blots, respectively. Bars represent smeared DNA without distinctive nucleosome laddering pattern. 



ethidium bromide-stained gels, suggesting that 72? RAP 1 
had an effect on chromatin structure at multiple loci in 
the genome. On the contrary, adding doxycycline to the 
control cells did not change the MNase digestion pattern 
(Supplementary Figure S8E and F). 

To examine chromatin structure at specific loci, we 
carried out Southern blotting using a VSG18 probe, a 
VSG2 probe, an mVSG397 probe and a TTAGGG 
repeat probe. In un-induced cells, most of the 
VSG18 DNA was not digested by a 5-min treatment of 
MNase, whereas in 72? RAP 1 -depleted cells, a 2. 5-min di- 
gestion had already degraded most of the VSG18 DNA 
(Figure 5 A, VSG18 hybridized blot, arrow heads). 
Depletion of Tb RAP 1 also had similar effects on the 
VSG2 (Supplementary Figure S8A), mVSG397 
(Supplementary Figure S8C), and the telomere chromatin 
(Figure 5B, black bars). Furthermore, we observed the 
same phenotype in PRi-C2 cells but not the control cells 
(Supplementary Figure S8B, E and F). The chromatin 
structure at the ES promotor region appears to be only 
mildly affected by Tb RAP 1 depletion (Supplementary 
Figure S8D), similar to what we observed in FAIRE 
analysis (Figure 4). 

We also performed the MNase digestion in BF cells. A 
normal laddering pattern was observed in ethidium 
bromide-stained gels after Tb RAP 1 was depleted, which 
was similar to that obtained before Tb RAP 1 depletion 
(Figure 5C and D; Supplementary Figure S9A and B). 
No significant difference was visible between un-induced 
and induced cells when the chromatin structure at the 



VSG18 (Supplementary Figure S9A), VSG2 (Figure 5C 
and Supplementary Figure S9B) and the telomere 
(Figure 5D) loci were examined. Control cells were not 
affected by doxycycline, either (Supplementary Figure 
S9C-E). Therefore, no significant change in chromatin 
structure was detected in BF cells on depletion of 
76RAP1. 

DISCUSSION 

TbRAPl as a key regulator of VSG silencing at different 
developmental stages 

Trypanosoma brucei has evolved sophisticated mechanisms 
to regulate VSG expression at different life cycle stages. At 
the insect stage, silencing both BF VSGs and mVSGs avoid 
their degradation by proteases in the midgut of tsetse (9). In 
BF cells, silencing of all m VSGs and all but one BES-linked 
VSGs is critical for successful antigenic variation and a 
long-term infection. Recent studies have unveiled several 
factors that affect BES-linked VSG silencing: Depletion of 
Tb RAP 1 in BF cells led to derepression of all BES-linked 
VSGs for typically several 10-fold (23); deletion of a histone 
methyl transferase, DOT lb, led to 10-fold higher of BES- 
linked VSG mRNAs in BF cells (20); depletion of 
7M)RC1, a component of ORC that is involved in DNA 
replication also led to several fold increase of BES-linked 
VSG mRNA in BF and PF cells (38). In addition, a number 
of factors including Imitation Switch (ISWI) (17,39), 
Nucleoplasmin-Like Protein (NLP) (40), Sptl6 (18) and 



7680 Nucleic Acids Research, 2013, Vol. 41, No. 16 



DAC3 (19) are necessary for proper regulation of BES 
promoter activities without affecting VSG silencing. It is 
clear that regulation of BES-linked VSG expression 
involves multiple mechanisms (14,15,17-20,23,36). On the 
contrary, regulation of mVSG expression appears to be 
entirely at the transcription initiation level (41,42), which 
is unusual in 7. brucei, where posttranscriptional regula- 
tion is predominant (43). However, mVSG silencing is 
poorly understood except that 7&ORC1, when depleted, 
led to mVSG derepression in PF cells (44). 

In this study, we found that depletion of 76RAP1 in PF 
cells has a similar effect as in BF cells — BES-linked VSG 
mRNA levels were increased several 10- to several 100- 
fold. In addition, 76 RAP 1 is also required for mVSG 
silencing at both BF and PF stages. Therefore, 76 RAP 1 
is required for silencing both BES-linked and mVSGs, the 
two types of VSG that can be expressed under normal 
physiological conditions, and 76 RAP 1 apparently has a 
much stronger effect on VSG silencing than other known 
VSG regulation factors. 

Regulation of VSG expression/silencing at BF and 
PF stages 

VSGs are expressed in a strictly monoallelic manner in BF 
cells and are silenced at the PF stage. In this study, using 
qRT-PCR analysis, we showed that the steady state 
mRNA levels of several BES-linked VSGs are much 
lower in PF cells than that in BF cells. A previous study 
indicated that the 3' untranslated region of the VSG gene 
renders its mRNA less stable in PF cells (36), which most 
likely is the reason for a lower VSG mRNA level in PF 
cells than in BF cells. 

Regardless of the steady-state levels of VSG mRNA in 
WT BF and PF cells, the changes in VSG mRNA levels 
induced by depletion of 76 RAP 1 reflect the silencing 
effect of 76RAP1. The strength of Tb RAP 1 -mediated 
silencing at the BES-linked VSG loci is significantly 
stronger in PF cells than in BF cells. Although we do 
not have all the proper cell lines to carry out a compre- 
hensive comparison for all BES-linked VSGs, comparison 
between our current and previous results suggest that 
most BES-linked VSGs are derepressed at a higher level 
in PF cells (several 10- to several 100-fold) than in BF cells 
(several to several 10-fold) (23). However, we noticed that 
different BES-linked VSGs were derepressed to different 
levels on depletion of 72? RAP 1. It is therefore possible that 
for some VSGs, the derepression levels are at more com- 
parable levels between BF and PF than other VSGs. 
Importantly, depletion of 72? RAP 1 affects chromatin 
structure significantly in PF cells, but not in BF cells, sug- 
gesting that modulating chromatin structure in PF cells is 
a mechanism of Tb RAP 1 -mediated VSG silencing. 

TbRAPl modulates chromatin structure in PF cells 

A number of recent studies have shown that modulation 
of chromatin structure is important for regulation of BES 
expression in BF cells. First, chromatin structures of the 
active and silent BESs are dramatically different (21,22). 
Second, several proteins with known functions in chroma- 
tin remodeling, including ISWI (17,39), DAC3 (19), Sptl6 



(18) and DOT lb (20) have been shown to play important 
roles in regulation of BES promoter activity. In this study, 
we showed that depletion of Tb RAP 1 led to a more 
loosely packed chromatin structure at all genomic loci 
tested and more significantly at subtelomeric VSG 
regions in PF cells, suggesting that TbRAPl -mediated 
VSG silencing involves modulation of chromatin 
structure. 

Depletion of TbRAPl in PF cells led to a striking 
change in the chromatin structure at subtelomeric VSG 
loci, including the minichromosome VSG671 locus, 
whereas the VSG671 mRNA level was not affected, 
indicating that the change in chromatin structure is not 
a consequence of elevated gene transcription, and that 
modulation of chromatin structure is a primary effect of 
76 RAP 1. VSG671 was not derepressed, presumably 
because it lacks a functional upstream promoter. 
However, for BES-linked VSGs and mVSGs, a loosened 
chromatin structure most likely is the cause or one of the 
reasons of their derepression. We noticed that changes in 
FAIRE values at some BES-linked VSG loci (including 
VSGs 2, 18, 15, 9 and 11) are slightly higher than that at 
the VSG671 locus (7 < 0.05). As all BES-linked VSGs 
were strongly derepressed on depletion of 76RAP1, it is 
possible that the elevated transcription of BES-linked 
VSGs may induce a further opening of the chromatin 
structure in addition to the initial effect resulting from 
the depletion of 76RAP1. 

At the insect stage, FAIRE and MNase digestion results 
indicated that Tb RAP 1 depletion affects chromatin struc- 
ture mildly at many loci throughout the genome, which is 
much broader than we previously anticipated. In both 
yeast and mammalian cells, RAP1 homologs have been 
found to locate at non-telomeric regions and act as tran- 
scription regulators (45-49), and yeast RAP1 has been 
found to help determine genome- wide chromatin structure 
(50). We previously observed that some 76 RAP 1 appears 
to localize at loci other than telomeres (23). Therefore, 
Tb RAP 1 may have a conserved function as its yeast and 
mammalian homologs in determination of genome-wide 
chromatin structure. 

To our surprise, we did not detect significant chromatin 
structure changes in BF cells on 76 RAP 1 depletion. Our 
observations suggest that the significantly more severe 
chromatin structure change and the significantly stronger 
gene derepression at BESs in PF cells than those in BF 
cells are linked. Currently, we favor several possible 
explanations (not necessarily mutually exclusive) for the 
different pheno types observed in BF and PF cells. First, 
the BES chromatin may have different structures in BF 
and PF cells, and Tb RAP 1 may be the responsible factor 
for this difference. An earlier work by Navarro et al. (51) 
found that T7 polymerase-mediated transcription from a 
chromosomally integrated T7 promoter is repressed along 
the entire length of the BES in PF cells, but not in BF cells. 
Therefore, the chromatin structure at the silent BES 
appears to be more tightly packed in PF cells. It is 
possible that 76 RAP 1 is the factor that maintains a 
more closed BES structure in PF cells. Second, as BF 
cells undergo antigenic variation and silent BESs need to 
be ready to switch to the active state, their silent state is 
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likely metastable, similar to what is often observed for 
telomeric silencing in yeast (24). Therefore, the BES may 
switch between the open and closed states more frequently 
in BF cells than in PF cells. In this case, 72? RAP 1 may still 
be involved in determination of chromatin structures in 
BF cells, but our techniques are not sensitive enough to 
detect a significant change owing to the fast switching 
status of the BES. Another possibility is that more 
factors are involved in modulating BES chromatin struc- 
ture and BES silencing independent of 72? RAP 1 at the BF 
stage. Therefore, depletion of 7&RAP1 in BF cells may 
still affect the chromatin structure, but a lesser degree of 
change or no quantifiable differences can be detected. 
Glover and Horn showed that an rRNA promoter- 
driven reporter gene targeted 5 kb from the telomere in a 
silent BES was expressed in PF cells, but not in BF cells 
(52), suggesting that additional factors are involved in 
silencing BESs in BF cells. A recent study suggested that 
depletion of histone HI led to a chromatin opening at 
multiple genomic loci only in BF cells, but not in PF 
cells (53), which is also consistent with this hypothesis. 

In this study, we further demonstrated that 72? RAP 1 is 
a key regulator of VSG expression in 7. brucei cells, and 
that VSG silencing appears to involve chromatin struc- 
ture. In PF cells, 7&RAP1 helps to determine the chroma- 
tin structure, particularly at the subtelomeric BES loci, 
which presumably contributes to Tb RAP 1 -mediated 
silencing. The effect of 76 RAP 1 on chromatin structure 
also appears to be more global than we anticipated. 

SUPPLEMENTARY DATA 
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